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ABSTRACT 

We measure the fraction of galaxies undergoing disk-disk major mergers (f^^^) at intermediate 
redshifts (0.35 < z < 0.85) by studying the asymmetry index A of galaxy images. Results are pro- 
vided for B- and ifg-band absolute magnitude selected samples from the Groth strip in the GOYA 
photometric survey. Three sources of systematic error are carefully addressed and quantified. The 
effects of the large errors in the photometric redshifts and asymmetry indices are corrected with maxi- 
mum likelihood techniques. Biases linked to the redshift degradation of the morphological information 
in the images are treated by measuring asymmetries on images artificially redshifted to a reference 
redshift of Zd = 0.75. Morphological K-corrections are further constrained by staying within redshifts 
where the images sample redward of 4000A. We find: (i) our data allow for a robust merger fraction 
to be provided for a single redshift bin, which we center at z = 0.6. (ii) Merger fractions at that 
z have lower values than previous determinations: = 0.045lQQit^ for A'Ib < —20 galaxies, and 

f^^P^ 0.03lt!]:J]j^^ for Mk, < -23.5 galaxies. And, (iii) failure to address the effects of the large 
observational errors leads to overestimating f™^^ by factors of 10%-60%. 

Combining our results with those on other i3-band selected samples, and parameterizing the merger 
fraction evolution as /™p''(z) = /^^''(O) (1 + ^)", we obtain that m = 2.9 ± 0.8, and f^'P^'iO) = 
0.012 ± 0.004. For an assumed merger time-scale between 0.35-0.6 Gyr, these values imply that only 
20%-35% of present day Mb < —20 galaxies have undergone a disk-disk major merger since z ^ 1 

Assuming a Kg-hand mass-to-light ratio not varying with luminosity, we infer that the merger rate 
of galaxies with stellar mass > 3.5 x 10^° Mq is SRm = l-Slo g ^ Mpc^^ Gyr"! at z = 0.6. 
When we compare with previous studies at similar redshifts, we lind that the merger rate decreases 
when mass increases. 

Subject headings: galaxies: evolution — galaxies:interactions — galaxies:statistics 



1. INTRODUCTION 

] Current A-cold dark matter simulations show that 
hierarchical halo mergers explain the build-up of dark 
' matt er structures in the Universe (jBlumenthal et aP 
; 119841: iSpringel et all 120051) . While the merger his- 
tory puts difficulties t o the formation o f disk 
■ galaxies (Ivan den Bosch' |200H lAbadi et all l2003t 
I lD'Onghia~ fc Burkcrt 2004), it is generally agreed 
that mergers are important in the formation o f mas - 
sive early- type galaxies (e.g., iDe Lucia et al.l l2006f ). 
[ Simulations suggest that gas- rich mergers can pro- 
I duce spheroidal systems Ce.g.. iNaab fc Burkertj l2003l : 
iBournaud eFall l2005t iHopkins et all l2008f) . ^i\e dis- 
sipationless spheroidal mergers ex plain more massive 
spheroids (e.g.. iGonzalez-Garcia fc van Albadal l2003l : 
iGonzalez-Garcia fc Balcellsl 120051: iNaab et 311120061 ). 
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Mergers rti a y also p lay a role in disk galaxy evolution: 
iLotz et all ()2008bl ) N-body simulations show that 
gas-rich major mer g ers can produce disk systems, while 
lEliche-Moral et al.l (|2006b[ ) show that minor mergers 
contribute to bulge growth in disk galaxies. Obser- 
vationall y, the size evolution of massiv e galaxies with 
redshift (|Truiillo et all l2007t Fvan Dokkum et al.l l2008f) 
and the luminosity density ey olutio ii of red sequence 
galaxies since z ~ 1 (jBeU et al.ll2"00l iFaber et al.ll2007f l 
rule out the passive evolution hypothesis and suggest 
galaxy mergers as an important process in galaxy 
formation. 

Despite their importance, the merger rate (Km, num- 
ber of mergers per comoving volume and time), the 
merger fraction (/m, fraction of mergers in a given 
sample), and their evolution with z, are poorly con- 
strained observationally. Merger fractions may be es- 
timated from statistics of close pairs, or from statis- 
tics of geometr ically-dis torted galaxies. Working with 
close pairs fe.g..lPatton e t al. 2000, 20"cil lLin et"al]|2004 
120081: be Propris et al.ll2005i. i2007. ) gives useful informa- 
tion on the progenitors of the merger, such as their mor- 
phologies, their mass ratio, or their relative color. On 
the other hand, these studies need spectroscopic sam- 
ples to avoid contamination by projection, which often 
leads to small sample sizes, although photometric red- 
shift samples a,re starting to be used for t his purpose 
(e.g., iBell et al." 2006': 'K artaltepe et al.ll206'7i: iRvan et all 
[20081: IHsieheraL,2008,l . 
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Methods based on morphological distortions, referred 
to in this paper as morphological merger fraction deter- 
minations, exploit the fact that, in the final stages of 
a disk-disk merger, the merger remnant is highly dis- 
torted, e.g., with high asymmetries , tidal tails, or double 
nuclei (IConselice et al.ll2003L l2008t ICassata et al.. ,2005;: 

m 



Kampczvk et all l2007t [Bridge et all 120071 : iLotz et all 
2008al hereafter LOS). The id entification of distorte d 



sources may be done by eye (e.g. JKampczvk et al.ll2007l ) 
or by automatic morph ological indices, such as the asym- 
metry index A ( Abraham et al.l 119961: IConselice et al.l 



I2000( ). and the G and M20 indices (|Lotz et al.l l200l ) 



When using automatic methods, the uncertainties in the 
morphological indices and in the photometric redshifts, 
coupled to the fact that distorted galaxies represent a 
small fraction of the total, lead to objects statistically 
'spilling over' from the most populated to the less pop- 
ulated bins, both in asymmetry an d in photometric red- 
shift. These effects were studied bv lLopez-Saniuan et al.l 
(|2008l hereafter LGB08), who used maximum likelihood 
(ML) techniques to quantify the errors, and to provide 
unbiased determinations of the merger fractions. These 
authors conclude that the use of classical, straight his- 
tograms to compute merger fraction evolution can easily 
lead to overestimating merger fractions by ~50%. Such 
errors lead to overestimating the importance of mergers 
in galactic evolution models, e.g., to reproduce galaxy 
number counts, or when comparing observational results 
with ACDM predictions. 

In this work we measure the morphological merger frac- 
tion at intermediate redshift (z ~ 0.6) using asym- 
metry index diagnostics. Previou s works find a merger 
fraction at that re dshift of ~ 0.07 (jConselice et al.ll2003l: 
Lotz et al.l ^ 2008a^■ although highe r (0.09. ICassata et all 



20051: 0.16. i Bridgc ct aU |200% . and lower (0.02, 
Kampczvk et a l. 2007) values have been reported. We 



apply the ML method developed in LGB08 to determine 
to what degree previous merger fraction determinations 
are affected by redshift and asymmetry measurement er- 
rors. 

The paper is structured as follows. In §|2]we summa- 
rize the GOYA data used in this paper, while in § [3] we 
describe the asymmetry index calculations and its vari- 
ation with redshift. In § |4j we review the ML method 
developed and tested in LBC08. The merger fraction 
values are presented in § [5l and we compare our results 
with those by other authors in § |6l Our conclusions 
are presented in § |7l We use Hq = 70 Km s~^ Mpc^^, 
fij\/ = 0.3, and fl^ = 0.7 throughout this paper. All 
magnitudes are in Vega system, unless noted otherwise. 

2. DATA 

We work with images and catalogs from the GOYA 
photometric survey, an imaging survey in preparation for 
the GO YA^ NIR spectroscopic survey with GTC/EMIR 
(jGuzm an 2003). For this paper, we focus on the Groth 
strip (GS) field, which is covered in six broadband fil- 
ters {U,B,V,I,J,Ks) over a common area of 155 square 
arcmin. The area covered is th at of the original Groth 
strip survey with HST/WFPC2 (iGroth at el.lll994l ). cen- 
tered on a = 14'^16"38"8 and S = 52°16'52" (J2000.0). 
IIST/WFPC2 imaging, which will provide the data for 



our asymmetry measure ments, have been extensively 
described elsewher e (e.g. iRatnatunga et al]ll995l 119991 : 
ISimard et al.ll2003 ). Exposure times for the images used 
here were 2800 s in F606W (Veoe) and 4400 s in F8UW 
(1314), and a typical depth of i^n — 25.0 is reached. 

j^s-band imaging is described in 

ICristobal-Hornillos et all ()2003f ). The Groth strip 
was covered with 11 pointings of the WHT/INGRID 
camera, with a pixel scale of 0.24 arcsec and seeing 
ranging from 0.6" to 1.1" FWHM. With typical exposure 
times of 5700 s, the median 3-cr depth is 20.5 mag. 
Imaging in the J band was similarly carried out with 
the WHT/INGRID camera. Exposure times were 1800s, 
leading to depths of 21.8 mag. 

U and i3-ba n d im aging are described in 
lEliche-Moral et al.l (|2006af l. The entire Groth strip 
was covered with a single pointing of the INT/WFC. 
Integration times were 14,400 s in C7 and 10,300 s in B, 
which led to ~3-cr depths of 24.8 mag in U and 25.5 
mag in B. The FWHM of the images was 1.3". 

2.1. The GOYA-GS Gatalog 

The parent catalog for the present study is a K^- 
selected catalog comprising 2450 sources, for which pho- 
tometry is provided in six bands, U through Kg. The 
limiting magnitude of the catalog, as derived from sim- 
ulations with synthetic sources, is Kg ^ 20.51 (50% de- 
tection efficiency). Multi-band photometry comes from 
images with matched PSF to 1.3" of FWHM. The cat- 
alog contains DEEP3^ spectroscopic redshifts (zspcc) for 
'-^ 600 sources and photometric redshifts (zp hot) for all. 
The latter were obtained with HyperZ (Bolz onella et al.l 
|2000( ). that also provided the 68% confidence interval 
for each Zphot- However, the probability distributions 
that describe the Zphot's are not symmetric, present 
non-analytic shapes and can have double peaks, while 
our methodology needs that errors in Zphot be Gaus- 
sian (§ m LGB08). In addition, we find that the con- 
fidence intervals given by HyperZ do not correlate with 
the differ ences between Zspec 's and Zphot's, a result also 
noted bv lOvaizu et al.l (|2008f ). Because of this, we used 
f^zphot = o-<5.(l + -Zphot) as Zphot error, where as, is the 
standard deviation in the distribution of the variable 
5z = (^phot - 2:spec)/(l + 2^phot), that is weU described 
by a Gaussian with Sz = —0.01 and as, ~ 0.07 (Fig. [T]). 
This procedure assigned the same error to sources with 
equal Zphot, but it is statistically representative of our 
sample and ensures the Gaussianity of Zphot errors in 
the merger fraction determination (§[!]), while 68% con- 
fidence intervals from HyperZ do not. 

HyperZ also yields the most probable SED of the 
source, that is used for degradation of the sources in the 
asymmetry calculation process (see §§ 13.21 and 13. 3p and 
for computing absolute magnitudes, necessary for doing 
the final selection of our sample (see the next section). 

2.2. Galaxy Samples 

We define two samples for our morphological analysis, 
selected in B- and Xg-band absolute magnitude, respec- 
tively. The i?-selected sample allows for comparisons 
with other studies in the literature, which often select 
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Fig. 1. — Histogram of the variable 5^ (see text for definition). 
The black solid line is the best fit of the histogram to a Gaussian 
with 1^ = -0.01 and cr^-^ = 0.07. 

their samples in that visual band. Comparisons are car- 
ried out in § 16.21 On the other hand, absolute Kg, mag- 
nitude is a good trace r of the stellar mas s of the galaxy 
(|Bell fc de Jond 120011 : iDrorv et al.ll2004D . This makes 
the galaxy selection less dependent on the instantaneous 
star formation, giving us a more nearly mass-selected de- 
termination of the merger fraction. 

The redshift range for our samples is determined as fol- 
lows. The highest redshift at which we can use asymme- 
try as a reliable morphological indicator without perform 
a morphological K-correction is Zup = 0.85 (see § I3.ip . 
In addition, because the ML method used in the merger 
fraction determination (see § H} takes into account the 
experimental errors, we must include in the samples not 
only the sources with zi < z^p, but also sources with 
Zi — 2ai < Zup in order to ensure completeness. Because 
of this, the maximum redshift in our samples, Zmax, must 
fulfill the condition Zmax ~ 2(T5^ (1 -I- Zmax) = 0.85, which 
yields Zmax ~ 1.15. We take as minimum redshift in our 
study Zmin = 0.2 because of the lack of sources at lower 
redshifts. This yields Zdown = •Zmin+2cr5^ (l+2:min) ^ 0.35 
to ensure completeness and good statistics. 

The limiting absolute magnitude of the samples is ob- 
tained by calculating the third quartile of the Mb and 
Mk^ redshift distributions in various redshift bins. We 
show those as filled circles in Figure [21 for Mb {upper 
panel) and M^^ (lower panel). The black solid curves 
are the least-squares fits to these points by a third de- 
gree polynomial. To be complete up to Zmax = 1.15 (in- 
dicated with vertical solid lines), the samples need to be 
restricted to Mb < -20 and Mk, < -23.5. These limits 
are indicated as horizontal dashed lines in both panels. 
The B — Ks color distribution of the two samples peaks 
at B — Ks — 3.6, consistent with the difference between 
both selection cuts, which is 3.5 magnitudes. This selec- 
tion yields 567 sources with Mb < —20 and < z < 1.15, 
and 505 with Mk, < -23.5 and < z < 1.15. We study 
the difference between both catalogs in § 15.11 

3. ASYMMETRY INDEX 

The automatic asyr nmetry index A is one of the CAS 
morphological indices (|Conselicell2003l hereafter C03) . It 
is defined as: 
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where Iq and Bq are the original galaxy and background 
images, /igo and Siso Sire the original galaxy and back- 
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Fig. 2. — (upper panel) and Mp^^ (lower panel) as function 
of redshift (gray dots). In both panels black dots mark the lim- 
iting magnitude at different redshifts. The black solid curves are 
the least-squares fits of the completeness points by a third degree 
polynomial. The horizontal black dashed lines are the magnitude 
cuts that we used to select our samples: Mb < —20 (upper) and 
Mk, < -23.5 (lower). 

ground images rotated by 180 degrees, and the sum spans 
all the pixels of the galaxy and background images. The 
index increases with the deviation of the galaxy image 
from point-symmetry. The background terms in equa- 
tion (III) are measured on a region of the frame free from 
known sources; it takes into account the level of asymme- 
try expected from the noise distribution in the sky pixels. 
F urther details on the a symmetry calculation are given 
in lConselice et al.l ()2000D . We use the A index to identify 
recent merger systems which are very distorted. On the 
basis of asymmetry measurements on images of nearby 
merger remnants, previous merger fraction determina- 
tions have taken a system to be a major merger remnant 
if its asymmetry index is A > Aj^, with = 0.35 (C03). 
Note that this criterion applies to disk-disk mergers only; 
spheroid-dominated mergers suffer much weaker morpho- 
logical distortions, hence they are missed by the asym- 
metry criterion just described. For high-redshift samples, 
the determination of A needs to be done on HST images 
to ensure high spatial resolution. In our case, we work 
with Vgoe and i^n bands. To increase the signal-to-noise 
we determined the asymmetry index of each source 
in the image Veoe + «8i4- 

3.1. Pass-band Restrictions to the Red.shift Range 

Galaxy m orphology depends on the band of obser- 
vatio n (e.g. iKuchinski et al.l 120001 : iTavlor-Mager et all 
|2007( ). In particular, when galaxies contain both old and 
young populations, morphologies may change very sig- 
nificantly at both sides of the Balmer/4000A break. The 
asymmetry index limit A^ — 0.35 was established in the 
rest- frame S-band (see C03). When dealing with galax- 
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Fig. 3. — Variation of fnpiz) with redshift. Black solid line is 
for B filter as reference, while gray dashed line is for U filter. The 



vertical black solid line indicates Zu 



■ 0.85. 



ies over a range of redshifts, in order to avoid system- 
atic passband biases with redshift, one needs to apply a 
so-called morphological K-correction by performing the 
asymmetry measurem ents in a band as clo se as possible 
to rest-frame B (e.g.. iCassata et al]|2005[ ). While cor- 
rections have been attempted for obtaining asymmetries 
in rest-frame B from asymmetry measurements in rest- 
frame U (Conselice et al. 2008), in the present study we 
stay within the redshift range where our images sam- 
ple rest-frame B. To determine the redshift ranges over 
which rest-frame B-band or [/-band dominates the flux 
in the observational Vqqq + igi^ filter, we define the func- 
tion 



/rf(z) 



/„°°Pv.(A/(l + z))PRF(A)dA 
/o°°PRF(A)dA 



(2) 



where Prf and Pvi are the transmission curves of the 
rest-frame reference filter and the Vgoe + *8i4 filter, re- 
spectively. In Figure[3] we show the function (black 
curve), and fuiz) {gray dashed curve). The redshift in 
which the [/-band starts to dominate the flux in the ob- 
served V(306 + *8i4 filter is z^p = 0.85 {vertical black solid 
line). We take this redshift as the upper limit for our 
study. 

3.2. Asymmetries at a Reference Redshift 

The asymmetry index measured on survey images 
varies systematically with the source redshift, due to the 
loss of spatial resolution and of source flux with z. The 
net result of such loss of information is a systematic de- 
crease with z of the measured asymmetry. Several pa- 
pers have attempted to quantify these effects by degrad- 
ing the image spatial resolution and flux to simulate the 
appearance a given galax y would hav e at different red - 
shift s in a given survey. IConselice et all (f2(M l200l : 
and ICassata et al] ()2005D degraded a few local galaxies 
to higher redshif ts, and found that ind eed asymmetries 
decrease with z. IConselice et al.l (|2003| ) also noted that 
this descent depends on image depth, an d that luminous 
galaxi es are less affected. In addition, IConselice et al.l 
()2005f l show that irregular galaxies (high asymmetry) 
are more affected than ellipticals (low asymmetry). A 
ze roth-order correc tion for such biases was implemented 
bv IConselice et al.l (2003, 2008) who apphcd a correction 
term AAz defined as the difference between the asym- 
metry of local galaxies measured in the original images 
and the asymmetry of the same galaxies in the images 



degraded to redshift z. Their final, corrected asymme- 
tries are A( = Aq + AA^, where Aq is the asymmetry 
measured in the original images. With these corrections, 
all the galaxies have their asymmetry referred to z = 0, 
and the local merger criterion A > Am = 0.35 is used to 
flag merger remnants. 

We improve on the above procedure by computing a 
correction term individually for each source in the cat- 
alog. Also, rather than attempting to recover z — 
values for A, we degrade each of the galaxy images to 
redshift Zd = 0.75; we then obtain our final asymme- 
try values A{ directly from the degraded images. With 
this procedure, we take into account that each galaxy is 
affected differently by the degradation, e.g., the asym- 
metry of a low luminosity irregular galaxy dramatically 
decreases with redshift, while a luminous elliptical is 
slightly affected. We choose Zd — 0.75 as our refer- 
ence redshift because a source at this (photometric) red- 



shift has Zd -t- (Tz 



0.85, that is, the probabil- 



ity that our galaxy belongs to the range of interest is 
~ 85%. Because we will work with asymmetries reduced 
to Zd = 0.75, we cannot apply the local merger criterion 
A > v4,n = 0.35. We redefine this criterion in ^ 13.41 

Only ~26% of the sources in the catalog have spectro- 
scopic redshifts, hence redshift information, coming pri- 
marily from photometric redshifts (§111), has large uncer- 
tainties. To account for the redshift uncertainty when de- 
riving the asymmetries at Zd = 0.75, we start from three 
different initial redshifts for each photometric source, 

^phot ~ ^phot — CTzphot I ^phot , and Zpj^i^^ = Zphot + i^Zphot I 

and degrade the image from these three redshifts to 
Zd = 0.75. Then, we perform a weighted average of the 
three asymmetry values, such that 

Af = 0.16^0.75 (Shot) + 0.16Ao.75(<hot) 

+ 0.68^0.75 (zphot), (3) 
where Ao.75{z) denotes the asymmetry measured in the 
image degraded from z to Zd = 0.75. When a spectro- 
scopic redshift is available, the final asymmetry is sim- 
ply A[ ~ ^o.75(^spoc)- We do not apply any degrada- 
tion to sources with z > 0.75, that is, we assume that 
^o.75(^ > 0.75) = ^0- Whenever a source is not de- 
tected after degradation, we remove it from the sample 
in spectroscopic redshift cases, and we do not use it in 
equation ([3]) in photometric redshift cases. 

To obtain the error of the asymmetry, denoted aAf , we 
average, in photometric redshift cases, the uncertainties 
of the three asymmetries following equation ([3]) , and add 
the result in quadrature to the rms of the three asymme- 
try values. The first term accounts for the signal-to-noise 
error in the asymmetry value, while the second term is 
only important when differences between the three asym- 
metry values can not be explained by the signal-to-noise 
first term. In the spectroscopic case we take as aAt the 
uncertainty of the asymmetry Ao.75(zspcc). 

The de gradation of the im ages was done with COS- 
MOSHIFT ( Balcells et al.ll2003l ). which performs repixela- 
tion, PSF change, flux decrease, and K-correction over 
the sky-subtracted source image. The K-correction for 
each source is computed through integrals of the best- 
flt SED obtained as output of HyperZ in the Zphot de- 
termination. The last COSMOSHIFT step is the addition 
of a random Poisson sky noise to the degraded source 
model. As a result of this last step, two COSMOSHIFT 
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degradations of the same source will yield different asym- 
metry determinations. We take the asymmetry of each 
degraded source, ^0.75(2), to be the median of asymme- 
try measurements on 10 independent degradations of the 
original source image from z \o z^ — 0.75. Note that, in 
photometric redshift cases, each final asymmetry value 
comes from three previous asymmetries (eq. [3]), so each 
Ai determination involves 30 asymmetry calculations. In 
all the cases the uncertainty in ^0.75(2) is the median of 
the 10 individual asymmetry errors. 

To check that the different final asymmetry determi- 
nations for sources with photometric and spectroscopic 
redshifts do not bias the asymmetry values, we com- 
pare the Ai of the 56 sources with Zspcc, Mb < —20, 
and 0.35 < Zgpcc < 0.75, denoted Af(spec), with the fi- 
nal asymmetries obtained from their corresponding Zphot 
and Equation [31 denoted Af(phot). The difference 
Af(spec) — At(phot) has an rms = 0.025, lower than 
the typical error in Af {'oaJ ~ 0.04). In fact, 90% of 
the sources have their >lf(phot) inside zLaAj{spec) and 
all sources inside ±2(TAj (spec). Therefore, we conclude 
that the A{ measured in sources with Zphot are equivalent 
to the A{ measured in sources with ^spec- 

The asymmetries A{ referred to Zd = 0.75 provide a ho- 
mogeneous asymmetry set that permits consistent mor- 
phological studies in the GS field. In § 16.11 we discuss 
the effects that the usage of degraded or non-degraded 
asymmetries has on the merger fraction determination. 

3.3. Asymmetry Trends with Redshift 

For a sample of galaxies over a range of redshifts, the 
statistical change of the measured asymmetries with z is 
the combined effect of loss of information (as shown in 
the previous section) and changes of the galaxy popula- 
tions. In contrast, the z evolution of Ai reflects changes 
in the galaxy population alone, given that morphological 
information in the images used to determine A{ is homo- 
geneous for the sample. We show here that the z trends 
of Aq and A{ are quite different. 

In the upper panel of Figure[4]we show the variation of 
Aq with redshift in the Mb < —20 sample, while in lower 
panel we see the variation of A{ for the same sample. In 
both panels, white squares are the median asymmetries 
in A(z) = 0.1 redshift bins, and the black solid line is the 
linear least-squares fit to the 0.5 < z < 0.8 points. In the 
Aq case the slope of the fit is negative, Aq oc — 0.015z, 
while in the A{ case the slope is positive, A{ cx 0.051z. 
In the first case, the negative slope reflects that the loss 
of information with redshift (negative effect on A) domi- 
nates over genuine population variations (positive effect, 
because at higher redshift galaxies are more asym met- 
ric, e.g.. ICassata et all 120051 : IConselice et al.|[2005() . In 
the second case, we have removed the loss of informa- 
tion term, so we only see population effects. We take 
as degradation rate, denoted Sa, the difference between 
both slopes, that yields AAb = 5^Az = — 0.066Az. In 
the ii's-limited sample we follow the same procedure and 
obtain AAk, = -0.060Az. 

3.4. Adjusting the Asymmetry-based Merger Criterion 

How do the trends shown in § 13.31 affect the merger 
criterion A^ = 0.35? This value was chosen in C03 to 
select the high asymmetry tail of the local asymmetry 
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Fig. 4. — Asymmetry vs. redshift in the Alg < —20 sample {gray 
dots in both panels). Upper, asymmetries of the sources measured 
on the original images. Lower: asymmetries of the sources mea- 
sured on images degraded to = 0.75. White squares in both 
panels are the mean asymmetries in 0.1 redshift bins. The black 
solid line is the linear fit to the mean asymmetries in the [0.5,0.8) 
redshift interval. 

distribution: most of the local galaxies with A > 0.35 
are merger systems. Note however that only ~ 50% of 
all local merger systems have A > 0.35 (C03, Fig. 9). 
If asymmetry systematically decreases with redshift, the 
local merger criterion needs to descrease as well to pick 
up the same distorted sources of the sample. This idea is 
supported bv lKampczvk et all ()2007f) . who visually com- 
pare distorted galaxies at z ~ 0.7 with their asymmetry 
values. They find that a third to a half of the galaxies 
with A > 0.20 ar e mergers, while 60% of merger systems 
have A > 0.20 (|Kampczvk et al.l 120071 their Fig. 5). 
Hence, although not all the sources with A > 0.20 are 
merger systems, the number of galaxies with A > 0.20 
are statically representative of the total number of merg- 
ers. We assume that the local merger criterion evolves 
with redshift as A^{z) = A„i(0) -I- Saz — 0.35 + 6az. 
With the two different values of Sa obtained previously, 
we infer that v4,n(0.75) ^ 0.30, and we use this merger 
criterion in the following. This variation is particular of 
our samples and HST/WFPC2 images: deeper images, 
and images with different angular resolution may show 
different trends. 

4. MAXIMUM-LIKELIHOOD MERGER FRACTION 
DETERMINATION 

Following iConselicd (|2OO60 . the merger fraction by 
morphological criteria is 

K N 

rmph ^ "^^^m 

iVtot + («-l)A^m' ^ ' 

where A^,„ is the number of the distorted sources in the 
sample with A > A,„, and iVtot is the total number of 
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sources in the sample. If k > 2 we obtain the galaxy 
merger fraction, f^^^, and k represents the average num- 
ber of galaxies that merged to produce one distorted rem- 
nant. If K = 1 we obtain the merger fraction, f^^^'- the 
number of merger events in the sample. We will use k — I 
throughout this paper. 

The steps that we follow to obtain the merger frac- 
tion are described in detail in LGB08. In this section we 
review the main steps. If we define a bidimensional his- 
togram in the redshift-asymmetry space and normalize 
this histogram to unity, we obtain a bidimensional prob- 
ability distribution defined by the probability of having 
one source in bin [zk,Zk+i) H [Ai,Ai+i), defined as pki, 
where index k spans the redshift bins of size Az, and the 
index I spans the asymmetry bins of size AA. We con- 
sider only two asymmetry bins split at Am, such that the 
probabilities pki describe highly distorted galaxies (i.e., 
merger systems) , while the probabilities pko describe nor- 
mal galaxies. With those definitions, the merger fraction 
in the redshift interval [zk, Zfc+i) becomes 



/'inph 
Jm,k 



Pkl 



PkO + Pkl 



(5) 



In LGB08 we developed a ML method that yields the 
most probable values of pki taking into account not only 
the z and A values, but also their experimental errors. 
The method is based on the minimization of the joint 
likelihood function, that in our case is 

L{zt,A,\p'kh'^z,,crA,) 



1- EE 

i ^ k I 

where 

ERF{r],i,k) = erf 



ERF(z,i,fc)ERF(A,i,0 



m - Vk+i 



erf 



Vi " Vk 



(6) 



(7) 



erf (x) is the error function, Zi and Ai are the redshift and 
asymmetry values of source i, respectively, <t^. and 
are the observational errors in redshift and asymmetry of 
source i, respectively, and the new variables pJcj = In(pfci) 
are chosen to avoid negative probabilities. Equation ^ 
was obtained by assuming that the real distribution of 
galaxies in the redshift-asymmetry space is described 
by a bidimensional distribution pki = exp(pj^.;) and that 
the experimental errors are Gaussian (see LGB08 for de- 
tails). Note that changing variables to p'^.^ = ln{pki), 
equation ([5]) becomes 



Jm,k ~ 



(8) 



LGB08 show, using synthetic catalogs, that the exper- 
imental errors tend to smooth an initial bidimensional 
distribution described by pki , due to spill-over of sources 
to neighboring bins. This leads to a ~ 10 — 30% over- 
estimate of the galaxy merger fraction in typical obser- 
vational cases. LOS find similar trends in their study of 
morphological merger fraction based on the M20 and G 
indices. LGBG8 additionally show that, thanks to the 
use of the ML method, we accurately recover the initial 
bidimensional distribution: the input and ML method 
merger fraction difference is ~ 1% even when experi- 
mental errors are similar to the bin size. That is, the 



ML results are not biased by the spill-over of sources to 
neighboring bins. 

We obtain the morphological merger fraction by ap- 
plying equation ^ using the probabilities pJcj recovered 
by the ML method. In addition, the ML method pro- 
vides an estimate of the 68% confidence intervals of the 
probabilities p'f.^, that we use to obtain the 68% 



confidence interval, denoted [ci^ 

f mph 



nph 1 Crtniph ] 



'm,fc 

- ,mp„, - „,„n ,- This interval 

is asymmetric because fc" is described by a log-normal 
distribution due to the calculation process (see LGB08 
for details). Note that, in LGB08, k = 2 is used in equa- 
tion ([?]), but the method is valid for any k value. 

4.1. Simulations with synthetic catalogs 

LGB08 show that the reliability of the ML method de- 
pends on factors such as the number of sources in the 
catalog, the mean experimental errors relative to the bin 
sizes, and the values of pki (bins with lower probabilities 
are more difficult to recover) . LGB08 also show that the 
probability distributions of pj,; must be Gaussian to en- 
sure the reliability of the method. We study the shape of 
probability distributions of p'^j by performing simulations 
with synthetic catalogs. 

We characterize the experimental catalogs with several 
parameters, which we use as input for the synthetic cat- 
alogs. These parameters are: the number of sources (n), 
the fraction of spectroscopic redshift sources in each red- 
shift bin (/fe,spoc), and the mean and the dispersion of Af 
errors in each redshift bin (cta ^, and CTo-^ Addition- 
ally, each experimental catalog have associated the p'j^^ 
probabilities that we obtain previously applying the ML 
method. To obtain these pj.; we fix the asymmetry bin 
size AA = 0.5, and only vary the redshift bin size Az, 
that is, Az is the only free parameter in this study. 

With the previous input parameters we create a syn- 
thetic catalog as follows: first we take n random sources 
distributed in redshift and asymmetry space following 
a bidimensional distribution defined by the probabilities 
Pkl = exp(pj.;). This process yields the input values of 
z and A, z-m^i and Ai^^i, of the n sources of our syn- 
thetic catalog, and the number of sources in each redshift 
bin, rifc. Next we apply the experimental redshift errors. 
For nkfk,spcc sources in each redshift bin we assume that 
the simulated redshift value is equal to the input value, 
Zsim.i = Zin.i, and assign it a constant standard devia- 
tion CTz^i^ i — 0.001. For the remaining nfc(l — fk,spec) 
sources in each redshift bin, the process is more compli- 



cated: we obtain the 



value as drawn for a Gaus- 



sian distribution with mean z^^^i and standard deviation 
(Jzsiry,,i = 0.07(1 -|- Zin^i). The process to obtain the simu- 
lated asymmetry values Asim,i is similar: these as drawn 
for a Gaussian distribution with mean Ain,i and stan- 
dard deviation ua^^^^ ^ ■ In this case, the value of o^A^i^-, i 
is a positive value also drawn for a Gaussian distribu- 
tion with mean cta and standard deviation cro-^ , so it 
depends on the redshift Zin,i of the source. 

In order to characterize the probability distributions 
that describe the p'^i, we generate a set of TV = 1000 
independent synthetic catalogs and apply the ML 
method to each catalog. We find that, due to the low 
number of highly asymmetric sources, the Gaussianity 
of the pj,; can only be ensured if we consider one redshift 



Major merger fraction at z — 0.6 



7 



bin at Zdown = 0.35 < z < z^p = 0.85 range (see § 12.21 
for details about these limits), that is, Az — 0.5. If 
we consider two redshift bins with Az = 0.25, the 
probability distributions of the p'f,i are non-Gaussian 
and we cannot ensure the reliability of the results. 



5. RESULTS 

On the basis of the arguments in §§[niand[31 we provide 
merger fractions for the redshift interval z € [0.35,0.85). 
In Table [T] we summarize the total number of sources, 
and the number of distorted sources {A{ > 0.30), with 
z G [0.35,0.85), for both classical counting and ML 
method. For the ML method, the number of sources is 
not an integer. Indeed, the ML method gives a statisti- 
cal estimation of the probability p^; = exp(pj,;) of finding 
one source in the redshift bin k and in the asymmetry 
bin Z, so the estimated number of galaxies in that bin, 
'T-fei.ML = n-totPkiAzAA (where ritot is the total number of 
galaxies in the sample), is not necessarily an integer. Ta- 
ble [T] shows that the number of distorted sources in the 
Xg-band sample is lower that in B-band sample. This 
result, which occurs for both classical and ML method 
determinations, is analyzed in more detail in § 15.11 In 
addition, the number of distorted sources given by the 
ML method is lower than that coming from the classical 
determination, but the total number of sources in the bin 
is higher. This is due to the fact that most of the sources 
in the samples (~ 70%) are in the range z £ [0.35,0.85), 
therefore more sources have spilled out of this bin due to 
redshift errors than viceversa. 

With the probabilities p'^,; and their confidence inter- 
vals given by the ML method, we obtain the following 
merger fractions for our B and Kg samples: 



f^^iz = 0.6, Mb < -20) = OMb^^ml (9) 

f^P\z - 0.6, Mk^ < -23.5) ^ Omit'^. fol. (10) 

In the next section we study the origin of the difference 
between both values, and we compare our results to other 
authors in § [6l 

5.1. Visible vs Near Infrared Merger Fractions 

We found that the merger fraction aX z — 0.6 in the 
Mk^ selected sample is a ~ 30% lower than in the Mb 
selected sample. Such trend had p r evious ly been noted 
in pair studies (jBundv et al]|2004t iRawat et aLH2008l ). 
To understand the origin of this difference, we study the 
nature of the galaxies in the range 0.35 < z < 0.85 that 
are not common to the two samples. For this discus- 
sion, we shall refer to galaxies only selected in B/K^ 
samples as the blue/red samples. The blue/red samples 
comprise 109/72 sources. As expected, the blue sam- 
ple comprises lower mass galaxies ((M^Jbiuo = —23.1 vs 
(M/fJrcd — —23.8), which are bluer than those from the 
red sample {{[B - ifsDbiuc = 2.6 vs {[B - K,]),^a = 4.4). 
We now show that the two samples have very distinct 
asymmetry distributions. In Figure [5] we plot the blue 
{open symbols) and red {filled symbols) samples in the 
A{-Mk^ plane. It is clear that asymmetries in the red 
sample are low, (v4„i)red = 0.10, with only 2 sources with 
A{ > 0.2 and none with A[ > 0.3 {black solid line). 
In contrast, the mean asymmetry of the blue sample is 
(^m)biuc = 0.15, with 11 sources with Af > 0.3 (-- 10% 




Fig. 5. — Asymmetry in function of Mk^ for blue (white dots) 
and red [black dots) samples (see the text for details). The black 
solid line shows the merger criterion Am = 0.30. 



of the blue sample). This result suggests that: (i) an im- 
portant fraction of the _B-band high asymmetry sources 
are low-mass disk-disk major merger systems that, due 
to merger-triggered st ar formation, have thei r B-band lu- 
minosity boosted up (|Bekki fc Shioval [20011) . enough to 
fulfill our selection cut Mb < —20. And, (ii) the objects 
not picked up in the B-selected sample are earlier types 
dominated by an spheroidal component which, when sub- 
ject to a major merger, do not distort sufficiently to be 
flagged as merger systems by our A-based criterion. 

6. DISCUSSION 
6.1. The Importance of ML Method 

Our merger fraction determination takes into account 
two potential biases. First, we artificially redshift the 
images to take into account the pixelation and signal-to- 
noise degradation with redshift (§§ [S21 ES])- Second, we 
use a ML method to take into account the experimen- 
tal errors in redshift and asymmetry (see LGB08). In 
Table [5] we compare merger fractions obtained with the 
same i3-band catalog but with different methods. First, 
we explore the values obtained if we simply count galax- 
ies without applying the ML method, i.e., if we mimic 
previous morphological merger determinations. In this 
case, tabulated uncertainties were estimated from Pois- 
son statistics. We see that, if we use non-degraded asym- 
metries Aq and the local merger criterion A > 0.35, 
we obtain a lower merger fraction than that from de- 
graded asymmetries A{ and the z — 0.75 merger criterion 
A > 0.30: 0.057 vs. 0.071 (20 vs. 25 merger systems, a 
~ 20% difference). Interestingly, when we apply the ML 
method over the two samples we obtain a similar merger 
fraction: 0.047 vs. 0.045, only a 4% difference. We com- 
ment on this fact later in this section. 

We also see that the effect of applying the ML method 
is more important when using A{ than when using raw 
asymmetries Aq: using ML with Af leads to a merger 
fraction decrease from 0.071 to 0.045, a ^ 60% difference, 
while, using Aq, it drops from 0.057 to 0.047, a - 20% 
difference. This difference is due to the fact that errors in 
Af are higher than errors in Aq: while Af values are su- 
perior tracers of the asymmetries of the sample galaxies, 
they have higher errors a a due to the uncertainties in the 
degradation process that transforms galaxy images from 
their original z to z^ = 0.75, and to the loss of signal-to- 
noise of those galaxies with z < z^. The mean Aq error in 
0.35 < z < 0.85 sources is WX^ — 0.02, while the mean Af 



8 Lopez-Sanjuan et al. 



TABLE 1 

Catalog Sources with z e [0.35, 0.85) 



Sample 


selection 


'^class 


riclass{^f > 0.30)'' riMh^ 


nML(Af > 0.30)^^ 


Mb < - 


-20 


352 


25 383.4 


17.4 


Mr, < 


-23.5 


313 


14 348.6 


10.8 



^ Number of merger systems in the bin.*^ Number of sources in the bin given 
by ML method. See §[5] for details. Number of merger systems in the bin 
given by ML method. See §[5]for details. 



TABLE 2 

Merger fractions in the Groth Strip at z = 0.6 



Method fS'^ {Mb < -20) /m'' (A/k, < -23.5) 





> 0.35=^ 




0.057 ±0.013 


0.035 ±0.011 


Ao 


> 0.35=^ H 


h ML'' 


047+°-''" 




At 


> O.SO'^ 




0.071 ± 0.015 


0.045 ±0.012 


Af 


> 0.30'^ 


f ML'' 







^ Using raw asymmetries Aq, and merger condition from C03 (z = 
0).^ Using ML method to determine the merger fraction.'^ Using 
z = 0.75 asymmetries A{, and merger condition at 2 = 0.75 13.31 1. 



error in 0.35 < z < 0.85 sources is WaJ = 0.05. As noted 
in LGB08, the bigger the experimental errors, the more 
the classical histogram-based methods overestimate the 
merger fraction. 

For the ifg-selected sample, the merger fractions de- 
pend on the usage of ML with similar trends to the B- 
selected sample. In the case we change from 0.035 
(11 merger systems) to 0.032, a ^ 10% difference, while 
in the Af case we change from 0.045 (14 merger systems) 
to 0.031, a ~ 45% difference. As in the B-band case, 
the ML method values are similar: 0.032 vs. 0.031, a 3% 
difference. 

The comparison of merger fractions in Table[5]indicates 
that, in the redshift range of our study, histogram-based 
merger fraction determinations may be overestimated by 
up to ^60%, and that such overestimates are readily cor- 
rected by using ML method. Once ML methods are used, 
the impact of working with asymmetry values normal- 
ized to a given reference redshift is less than 4%. Note 
that the effect of the experimental errors depends on the 
sample: our near infrared merger fraction is less affected 
by observational errors, but we do not know if this is a 
general trend, or a peculiarity of our samples. 

In summary, if we do not use ML method to take into 
account the effect of observational errors, we overesti- 
mate the merger fraction by 10%-60%, in good agree- 
ment with the expected ^10%-30% obtained in LGB08 
with synthetic catalogs. On the other hand, as long 
as ML method is used, the effect of the pixelation and 
signal-to-noise degradation with redshift is less impor- 
tant, only a ~ 4% effect. This fact, however, depends 
on the sample. A similar study in GOODS-S reveals 
that loss of information with redshift is also an impor- 
tant bias that needs to be treated to ensure accurate 
results (Lopez-Sanjuan et al., in prep.). 

6.2. Comparison with Previous Morphological B-hand 

Studies 



In this section we compare our morphology-based 
merger fraction to other determinations in the litera- 
ture. We restrict the comparison to works with B- 
band selected samples with well established luminosity 
limits, given the dependence of the merger fraction on 
the selection band, § 15.11 and given that more luminous 
- B-band samples tend to have higher merger fractions 
(jConselice et al.l[2003[ ). We exclude studies ba sed on pair 
statis tics, due to the progenitor bias (L08; iBell et aP 
|2006( ): each distorted galaxy in our sample is the final 
stage of the merger of two less luminous/massive galax- 
ies. For example, assuming a 1:1 merger and neglecting 
star-formation, we need merge two Mb = —19.25 galax- 
ies to obtain one high asymmetric Mb — —20 source. In 
addition, highly merger-triggered star-formation affects 
the S-band luminosity of the distorted s ource from -1-0.5 
(dust y merger) to —1.5 magnitudes (see lBekki fc Shioval 
[2OOII for details). 

In the following sections, we pay attention to how 
the information degradation and the experimental er- 
rors have been treated in previous works. The effect of 
the information degradation was addressed in all pre- 
vious morphological studie s at i n termediate redshifts 
(e.g.. IConselice et alj [20031 I2Q08I: ICassata et all 120051: 
iKampczvk et al.ll2007t L08), while the effect of experi- 
mental errors was studied in detail only by LOS. Using 
synthetic catalogs in the same way as in this paper, 
and in LGB08, they found similar trends: the experi- 
mental errors tend to overestimate the merger fraction. 
However, LOS did not apply any correction for this ef- 
fect. Our work is the first in which this important bias 
is corrected. 

6.2.1. The Merger Fraction at z ^ 0.6 

The mo st direct comparison is the work of 
iConselice^ t al. ( 20 0l) , wh o use t he same asymmetry in- 
dex as us. [Conselice et al.l ()2003[ ) lists a merger fraction 
/mph(2 ^ jVfs < -20) = 0.07. Their value is re- 
ally an upper limit coming from 1 merger detection in 15 
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galaxies, i.e., /™p'' < 0.16 assuming a Poisson distribu- 
tion. 

Other asymmetry works at intermediate redshift 
have been used di fferent selection criteria than us: 
ICassata et al.l (j2005) obtain a merger fraction /™p'^(z = 
0.75) = 0.088tS];[!26 in a tuk, < 20 selected sample. 
This value is a factor two higher than ours; we sus- 
pect that properly accounting for the experimental errors 
would reduce this discrepancy. Bridge ct al. (2007) per- 
form their asymmetry study on a 24um-selected sample 
{LiR > 5.0 X 101° Lq), finding /™p'X^ = 0.75) - 0.16. 
Such high value are expected from the fact that strong 
star formation occ urs in morphologically-distorted galax- 
ies (jSanders et al.| [l988). 

Two other studies give the morphological merger frac- 
tion using methodolog i es rel ated but not identical to 
ours: iKampczvk et al.f ()2007[ ) measured the fraction of 
visually disturbed galaxies at z '--^ 0.7, finding = 
0.7, Ms < -19.15) = 0.024, while LOB used the morpho- 
logical indices G and M20 to obtain /™p''(2; = 0.7, Mb < 
-19.75) = O.Oltom- In next section we use LOS data to 
constraint the merger fraction evolution, showing the dif- 
ficulties entrained when combining merger fractions from 
different methodologies. 

In summary, previous morphological works in B-band 
selected samples suggest that /™p'^(z ~ 0.7) ^ 0.07, 
higher than our value of 0.045. We recover their value 
(/™P^ = 0.071) when mimicking their methodology, i.e., 
when we apply the merger condition Af > 0.30 with- 
out using ML method. This suggests that merger frac- 
tion values around f^^^^iz - 0.7) ~ 0.07 may be over- 
estimated by ~50% due to not properly accounting for 
sources spilling over to neighboring bins. 

6.3. Morphological Merger Fraction Evolution 

We now constrain the evolution of the merger fraction 
with redshift by combining our result at z ^ 0.6 w ith 
those of L08 at 2 = 0.9, and iDe Propris et all (|2007f) at 
z = 0.07. 

To combine our results with LOB we check, first, 
that their luminosity-dependent sample selection, Mb < 
— 18.83 — 1.3z'^, matches our luminosity selection at z ~ 
0.9. We then note that LOB use different morphological 
indices than ours, namely G and M20. This introduces 
a time factor AT„i that relates the time that distorted 
sources fulfill the merger criteria in each of the methods: 
ATm — Tni.A/Tni,GM2o- The different merger time-scales 
can be constrained with merger simulations. iLotz et "all 
(|2008bf ) performed N-body/hydro-dynamical simulations 
of equal-mass gas-rich disk mergers and studied the time- 
scales Tai^A and Tm,GAf2o as a function of the merger pa- 
rameters. Their results suggest that ATm = 1.5 ±0.5. In 
addition, and taking into account the results of § 16. li we 
estimate a 20% overestimate on the merger fraction due 
to morphological index errors ((ctg) ^ 0.02). With all 
these considerations, the wet merger fraction at z = 0.9 
from LOB, f^^^iz = 0.9) = 0.06 ± O.OL* becomes 



mph,L08 



(z = 0.9, Mb < -20) = 0.075: 



-0.03 
-0.03 



(11) 



^ This luminosity cut selects Lg > OALg galaxies and takes 
into account the evolution of Lg with redshift l|Faber et al.ll200^ . 

* LOS find 41/685 blue (wet) mergers and 15/685 red (dry) merg- 
ers at z = 0.9. 




Fig. 6. — Morphological merger fr action in function of red shift 
for Mg < —20 galaxies. The data are lDe Propris et al.l 1 120071 . gray 
triangle), LOS {gray dot), and this work {white dot). The black solid 
line is the best fit to data, /m''''(2, Mg < -20) = 0.012(1 + z)2-9. 



in our methodology. The error is dominated by the un- 
c ertainty in the time-scal e factor AT^. 

iDe Propris et al.l ()2007f) provide a suitable low-redshift 
determination of the merger fraction using asymme- 
tries, but their result cannot be directly combined with 
ours since they exte nd to fainter absolute magnitudes, 
namely Mb ^ — 19. IConselice et al.l ()2003[ ) results sug- 
gest that f^^'^iMB < -20) ~ 1.5/^p'^(Mb < -19) by 
asymmetries, while pair statistics at that redshift also 
suggest that f^^HMfi < -20) > /iT''(Mb < -19) 
(jPatton fc Atfieldll2008l)'^. Beca use of this, we apply the 
C03 factor to iDe Propris et al.l ()2007[ ) results, and esti- 
mate that 



^mph,P07(^ ^ 0.07, Mb < -20) - 0.014™_ 



(12) 



We parameterize the merger fraction as /™p'^(z) — 
f^P^iO) (1 + z)" (e.g. JLe Fevre et al .1120001) . and per- 
form a weighted least-squares fit to the data, obtaining 
that 

f^P^\z,MB < -20) = (0.012±0.004)(l-f-z)2-3±° ^ (13) 

Figure [6l shows the three merger fractions from 
IDe Propris et al.) (j2007l . gray triangle), LOB {filled cir- 
cle) and this work {open circle). The black solid 
line is the weighted least-squares fit to the data, 
equation ([13]). This result supports the idea of im- 
portant evolution of the merger fract ion (m > 2) 
from z ^ 1 to the presen t (e.g., |Le Fevr e et al. | 
2000| iLaverv e t al. 2004; Cassa ta et al.|[2b 05: Be ll et al.l 
2006[ IKam pczvk ct al. 120071: iKartaltc pe et all l2007r 
Rawat et al. 2008; Hsieh et al.l [2008.1 . This evolution. 



however, do not imply that disk-disk major mergers 
were important in galaxy evolution since z 1: the lo- 
cal merger fraction is low, ^ 0.01, and, despite of the 
m ~ 2.9 evolution, the merger fractions remain below 
0.1 up to z ~ 1. It is only for z > 1 that our high expo- 
nent translates into a high merger fraction: f™^^ ^ 0.2 
at z ^ 1.5, extrapolating our fit. At those redshifts, ma- 
jor disk-disk mergers are sufficiently fre quent that they 
may be importa nt for galaxy evolution (Conselice '2006|; 
IConselice et al.l i2008). We study the importance of this 
type of mergers since z = 1 in the next section. 



5 IPatton fc Atfieldl 1I2OO8I ') perform their study in Mr (SPSS) 
absol ute magnitude samples. We take B-r = 1.25 IFukugita et al.l 
119951) to obtain the equivalent Mg samples. 
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It is worth recalling that our m — 2.9 do not necessary 
imply an important decrease with cosmic time of the total 
merger fraction since z = 1. Our merger criterion is only 
sensitive to disk- disk major merg ers. The pair major 
merger studies of iLin et al.1 HQQI); and Ide Ravel et all 
HqOSi) show that the wet merger fraction evolution is 
higher than the total (dry-|-wet) merger fraction evolu- 
tion. Our study cannot rule out an increasing importance 
of dry mergers at z < 1. 

6.3.1. Major merger remnants since z ^ 1 

Following iPatton et al.l (|2000l ). the fraction of present- 
day Mb < —20 galaxies that have undergone a disk-disk 
major merger since a given redshift is: 



/rc 



^ 1 



- 0.5/-P^z,) 



(14) 



where fai^^{z) is the merger fraction at redshift z, pa- 
rameterized as f^^^{0) (1 -I- z)™, Zj corresponds to a 
look-back time of t = jTm.yi, and N is the number of t 
steps since a given redshift z. We use the merger frac- 
tion parameters obtained in the previous section, and 
take two values for the merger time-scale: T^^a — 0.6 
Gyr (iLotz et al.l l2008bl . from N-body/hydro-dynamical 
equal-mass merg er simulations), and T,„,a = 0.35 Gyr 
([Consehcg l2006l from N-body major merger simula- 
tions). 

We find that the disk-disk major merger remnant frac- 
tion since z = 1 is /rem 20% for Tm a = 0.6 Gyr, and 
/rem ~ 35% for Tm,A = 0.35 Gyr. 

6.4. The Major Merger Rate at z = 0.6 

We now use the derive merger fraction for the Kg sam- 
ple to study the merger rate and its variation with galaxy 
mass. We define the major merger rate '^ra{z^ MkJ as: 



3?, 



,(z, MkJ - n(z, M^J/r^^, MkJT-^a, 



(15) 



where n{z,MK^) is the comoving number density of 
galaxies at redshift z brighter that Mk^, and T^^a is 
the merger time-scale, which we take as Tm,A = 0.35 — 
0.6 Gyr (^ I5XTD. T o obtain n{z,MK^ we use the 
ICirasuolo et al.l ( 2008[ ) luminosity function, that yields 
n(0.6,- 
obtain 



23.5) = 0.0025 Mpc"-^. With these values, we 



3fim(0.6,-23.5) 



l.Qtr, X 10- 



Mpc Gyr 



(16) 



The error takes into account both merger fraction and 
merger time-scale uncertainties. 

Adopting a constant mass-to-light ratio in the Mk. 
band of M^/Lk = 0.7 (iDrorv et al.ll2004l : lArnouts et all 
[200l . our Mk^ < —23.5 luminosity selection corre- 
sponds to a Mi, > 3.5 X 10^° Mq mass selection. Our 
inferred merger rate is shown against galaxy stellar mass 
in Figure[71 and listed in TablelSl together with literature 
values for the major merger rate at z = 0.6 at various 
limiting masses. Errors take into account both merger 
fraction and merger time-scale uncertainties. We can see 
that the merger rate at z = 0.6 decreases with mass. We 
find that the variation with M^, is well described by 



TABLE 3 

Merger Rates at z 


= 0.6 


Reference 


Stellar mass 

(Me) 


{10-'^ Mpc-3 Gyr-i) 


Conselice et al. (2008) 
Conselice et al. (2008) 
Conselice et al. (2008~) 


108 
109 

lOio 


46.9 ± 23.7 
26.4 ±14.7 
4.6 ±4.6 


This work 

Bell et al. (2006) 


3.5 X 10i° 
6 X lO^"^ 


^•"-0.6 

1.1 ±0.2 




I I I 

9.0 9.5 10.0 10.5 11.0 
log (M./Mq) 

Fig. 7. — Major m erger rate in function o f stellar mass a.t z = 0.6. 
The d ata are from [Conselice et al.l I I2008I . gray dots), IBell et all 
l|200a . gray square, error bars are smaller than the symbol si ze) , 
and this work {white dot). The solid line is the best equation (117( 1 
fit to the data. 

where M = log(Af*/Mo), and Mq, ^q, and /3 are param- 
eters to fit. The best fit to the data yields Mq = 
3.2 X 10^ Mq, % = 5.3 X 10- 3 Mp c~^ Gyr-\ and f3 = 
-0.36. Note that iBell et al.l (|2006f) merger criterion is 
sensitive to both disk and spheroid mergers, w hile asym- 
metry studies are only to disk-disk mer gers (IConselicd 
I2006f ). This implies that the data from B ell et all ()2006[ ) 
must be higher than the disk-disk merger rate. If we re- 
peat our analysis without the 6x 10^ Mq point, we obtain 
Mo 6.3 X IQi" Mq, % = 4.8 x lO^^ Mpc"^ Gyr"\ 
and P — —0.45. With these values, the inferred disk- 
disk m erger rate at 6 x 10^° Mq is 80% of the Bell et al.l 
(j2006f ) value. 

In addition, our infered major merger fraction for 



+0.013 
-0.009' 



3fim(M) = 3fioe'3^^' 



(17) 



M^ > 3.5 X 10^° Mq galaxies, /^p*^ = 0.031 
agreement with the visual estimate of the morphological 
major merger fraction of A/* > 2.5 x 10^" Mq galaxies in 
the r ange z e [0.34,0.8], f^^^ - 0.02 ±0.01 (|Jogee et al.l 
[200l . 

7. CONCLUSIONS 

From Groth strip survey images, we provide a robust 
determination of the disk-disk major merger fraction 
based on morphological criteria. We have quantified and 
corrected for the bias due to varying spatial resolution 
and image depth with redshift, by artificially redshift- 
ing the galaxy images to a common reference redshift of 
Zd — 0.75. More importantly, we successfully accounted 
for spill-over of sources to neighboring bins caused by the 
errors in asymmetry indices and in Zphot, through the 
use of a ML method developed in LGB08. The merger 
fractions for the i?-selected and Xg-selected samples are, 
respectively. 



0.6, Mb < -20) 



n45+° °i'' 

"•L'^3-0.011' 



Major merger fraction a,t z = 0.6 
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fZ^'iz = 0.6, Mk, < -23.5) = 0.03i;°;°J3. 

The effect of the experimental errors is the dominant 

observational bias in our study: without the ML method 
we overestimate the galaxy merger fraction by up to 60%. 
In comparison, the loss of information with redshift only 
biases the results by ~ 4%, as long as ML is used to 
account for the experimental errors. 

Parameterizing the merger fraction as f^^^ = 
/™P'^(0)(1 + 2)'", we obtain m = 2.9 ± 0.8, /,'^p''(0) = 

0. 012 ± 0.004. With these values, we infer that only 
20 — 35% of present-day Mb < —20 galaxies have un- 
dergone a disk-disk major merger since ~ 1. 

We use the M^^-band merger fraction to obtain the 
major merger rate at z = 0.6. Assuming a constant 
mass-to-light ratio, we obtain 3?i„(0.6, 3.5 x 10^^ Mq) = 

1. e^o g X 10"^ Mpc"^ Gyr"^ We compare our results 



with previous merger rates at that redshift, showing that 
the merger rate rapidly decreases with mass, such that 
the rate at M^, = 10^°-^ Mq is 10 times lower than that 
at = 10^ Mq. 

We dedicate this paper to the memory of our six 
lAC colleagues and friends who met with a fatal acci- 
dent in Piedra dc los Cochinos, Tencrifc, in February 
2007, with a special thanks to Maurizio PannicUo, whose 
teachings of python were so important for this paper. 
This work was supported by the Spanish Programa Na- 
tional de Astronomia y Astroffsica through project num- 
ber AYA2006-12955. 

Facilities: HST (WFPC2); ING:Herschel (INGRID); 
ING:Newton (WFC) 
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